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ABSTRACT 

A  model  is  described  which  incorporates  the  generation  and  ionization 
of  atoms  during  the  vaporization  of  individual  solute  particles  in  flames. 
Particles  initially  larger  than  a  critical  radius  have  two  vaporization 
rate  constants  because  of  a  change  in  the  vaporization  mechanism  related 
to  size.  Consequently,  the  mathematical  formulation  of  the  model  consists 
of  a  series  of  three  equations  which  are  used  to  curve-fit  different  portions 
of  emission  intensity-time  profiles  in  an  air-acetylene  flame.  Values  for 
the  vaporization  rate  constants  for  three  alkali  chlorides  found 
for  this  fitting  process  are:  Na  (244  vim2/s,  127  pm/s),  K  (650  um2/s, 

:q>y  nm/s),  and  Cs  (585  pm?/s,  847  pm/s) .  First -order  ionization  rates  arc 
58  s'1  for  Na,  268  s'5  for  K,  and  1128  s"1  for  Cs. 
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In  flame  spectrometry,  a  wet  aerosol  of  sample  solution  must 
undergo  dcsolvation  and  vaporization  to  form  the  free  atoms  needed  for  analysis. 
During  desolvation,  the  solvent  associated  with  the  sample  evaporates,  leaving 
behind  the  dry  solute  particle.  The  desolvation  process  has  been  studied 
previously  [1],  Vaporization  begins  immediately  after  desolvation  of  the  par¬ 
ticle  is  complete:  apparently  both  atomic  and  molecular  species  are  volati¬ 
lized  from  the  surface  of  the  dry  particle.  Vaporization  in  flames  has  been 
reported  [2]  to  occur  through  either  a  heat  or  mass-transfer  controlled  process, 
depending  on  the  composition  of  the  solute  particle.  In  this  communication 
a  model  is  proposed  which  closely  fits  experimental  data  on  atom  formation 
gathered  in  our  laboratory  using  the  isolated  droplet  sample  introduction 
method  [3].  Simplified  calculations  based  on  this  model  provide  approximate 
but  useful  values  for  several  rate  constants  important  in  the  particle 
vaporization  and  atom  ionization  processes. 

Tlll-'ORY 

'WWW 

The  conclusion  reached  in  previous  studies  of  vaporization  [2]  was 
that  there  arc  two  equations  wh ich  describe  the  time-dependent  radius  of  a 
solute  particle  vaporizing  in  an  analytical  flame.  Mien  the  particle  di¬ 
ameter  is  much  larger  than  the  mean  free  path  of  the  free  atoms  (4]  (iafge 
particle  model),  the  equation  for  the  radius,  r,  of  the  dry  particle  as 
a  function  of  time,  t,  is  similar  to  that  for  desolvation  [1,2]: 

-  kVit)’/-'  t  rV/ky, 


r  =  (r„7 


0) 


s 


where  ro  is  the  initial  particle  radius  and  k^  is  the  vaporization  rate 
constant  for  the  large  particle  model.  For  particles  much  smaller  than 
the  mean  free  path  of  the  released  atoms  [4]  (small  particle  model),  the 
equation  becomes: 

r  "  r«  -  kV2t  t  <  r»/ky2  (2) 

where  ky2  is  the  vaporization  rate  constant  for  the  small  particle  model. 

Both  equations  are  expected  to  he  valid  at  particular  periods  during  the 
life  of  a  single  solute  particle  except  when  the  initial  diameter  of  the 
particle  is  much  smaller  than  the  mean  free  path  of  the  liberated  atom.  In  the 
present  case,  the  mean  free  path  is  of  the  order  of  0.3-1. 5  pm,  as  discussed 
later.  In  this  latter  case,  only  the  small  particle  model  would  be  valid. 

The  number  of  atoms,  F,  which  are  vaporized  will  be  proportional  to  the 
mass  vaporized  (m)  which,  in  turn,  is  proportional  to  the  cube  of  the  particle 
'radius  (2]: 


Nm  _  4jnpNc  f  3  3) 

l  3M  '■  o  r  > 


(3) 


where  N  is  Avogadio’s  number,  M  is  the  atomic  weight  of  the  analyte,  e  is  the 
fraction  of  mass  entering  the  vapor  phase  as  free  atoms  (atomization  efficiency), 
and  p  is  the  density  of  the  analyte  atoms  in  the  condensed  pliasc.  Combining 
equations  (1)  and  (3)  for  the  large  particle  model  yields 


and 


for 


..  4iTpNc  r_  3 
*•  “  ~ 3M  lr« 

the  small  panicle 

i;  , 

M  o 


-  •  kv,t>3/!> 

model  equations  (2) 


0  :  t  <  r0*/kvl 
and  (3)  give 

0  5  t  <  r0/kw 


CO 


(5) 
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In  the  vaporization  of  many  solutes,  particularly  those  of  the  alkali 
metals,  atom  loss  through  ionization  must  also  be  considered.  For  these 
elements,  a  situation  exists  in  which  there  is  the  formation  of  the  atomic 
species  through  vaporization  of  the  particle  and  the  simultaneous  loss  of  the 
species  by  ionization.  Previous  studies  of  alkali  metals  [5,6]  have  found  the 
ionization  rate  to  he  first  order  with  respect  to  the  metal.  Of  course,  an 
observed  loss  of  alkali  atoms  through  ionization  would  not  follow  strictly 
first-order  kineiics  because  of  the  simultaneous  occurrence  of  ion-electron 
recombination,  an  event  which  opposes  the  ionization  process.  However,  at 
short  times  (i.e.,  within  two  ionization  time  constants,  2k^  *)  after  the 
completion  of  alkali  particle  vaporization,  the  contribution  of  ion-electron 
recombination  to  the  ionizxtiou  process  appears  to  be  small  and  constant, 
being  caused  principally  by  intrinsic  flame  electrons.  During  this  period, 
the  overall  vaporization-ionization  process  for  the  large  particle  model  can 
.  be  represented  by  t he  following  differential  equation: 


dF  =  VrrjaNf 
dt  '  3M 


2  kv,(ro2 


k]F 


(6) 


where  F  is  the  number  of  atomic  species  at  any  time  (t)  and  k^  is  the  ionization 
rate.  To  solve  tins  equation,  it  can  be  rewritten  in  the  following  form: 


rck(t  -_WNc  f  3^,^  .  kvit)lA  ek,tdt  ,  c 

If  the  expoiu i* t  i a  1  K  nn  in  tin-  integral  is  expanded  in  a  series  and  only  the 
fir*  I  term  i  ■  nl  i"ned,  l  lie  :*pi  me  iitc  solution  of  equation  (?)  is: 


i 


F  =-$r  Cro3  -  (r02  -  ku,t)^]e'kit 


The  number  of  free  atoms  (F)  can  be  related  to  the  radiant  intensity  (I) 


they  emit  by 


I  =  F(othvA)  (4tt)  “ 1  (W  sr-1) 


where  u  is  the  fraction  of  atoms  which  are  excited  (Boltzmann  ratio),  and 
hv  the  energy  and  A  the  emission  rate  (Einstein  coefficient)  for  the  observed 
transition.  Therefore,  eq.  8  can  be  expressed  as 


I  =  K[r03  -  (r02  -  kVlt) ^fe]e“kit  0  <  t  <  r02/kVl 


(10a) 


I  «  Krn 3 


t  >  r02/kVi 


(10b) 


where  K  combines  the  pro -parenthetical  constants  in  cq.  8  and  the  free-atom- 


to- intensity  conversion  of  eq.  9. 


4ttpMc 
K  =  3M 


ic  /«hvA\ 
\4tt  )  = 


pNemhvA 

3M 


In  actual  measurements,  the  radiant  intensity  (I)  is  seldom  measured; 
instead,  relative  emission  intensity  values  arc  recorded.  Consequently, 
the  experimental  value  for  K  which  is  obtained  includes  a  number  of  scaling 
factors  (including  measurement,  solid  angle,  optical  transmission,  detection 
cl  tic iency ,  electronic  gain,  etc.)  and  can  best  be  derived  empirically. 

'the  empirical  value  will  In'  <h  •  ijai.iM-d  K’  in  the  later  treatment. 

A  similar  mathematical  ti'.ituriit  of  the  small  particle  model  is  pos- 


dF  _  *nrpNe  r^.  ,  ,  2,  , 

dt  "  3M  l 3 k Vz  ( r o  -  kV2t)2]-k^F 


Fe  1 1  = 


t  _  ^pNe  / 
-~3M- J  3k 


3kV2(r0  -  kV2t)2  dt  +  C 


p  -  ~-N-r- [ ru 3  -  (r0  -  kV2t)3]e"kit 


I  -  K[r03  -  (r0  -  kVit) 3]e"ki1 


0  <  t  <  r0/k. 


o 3e~ki^  t  >  r0/kv 


I  =  Kr 


In  the  most  general  case,  the  vaporization  and  ionization  processes 
can  be  combined  for  both  large  and  small  particles  in  a  series  of  equations 
which  relate  intensity  to  the  initial  particle  radius  and  to  a  parameter 
we  have  termed  the  critical  radius,  r  .  The  critical  radius  is  that  hypo¬ 
thetical  particle  size  at  which  control  of  vaporization  in  the  model  appears 
to  change  from  the  large  particle  to  the  small  particle  mechanism.  The 
relationship  of  the  critical  radius  to  the  critical  time,  tc»  and  the 
time  of  complete  vaporization,  t can  be  found  from  equations  (17)  and  (18). 
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tc  =  (r02  -  rc2)/kVl  r0  >  rc 


(17) 


tf  =  (r02  -  rc2}/kv,  +  rc/kV2-  r0  >  rc 


(18) 


The  three  equations  required  to  describe  the  entiie  vaporization-ioniza¬ 
tion  process  for  the  combined  model  are: 

I  =  K[ro3  -  (r02  -  kVlt) 0  <  t  <  tc  (19) 

1  =  K[r03  '  (rc  “  kv2  tc  -  t  *  tf  (20) 


I  =  Kr03  e"1'^  t  >  tf 


(21) 


Equation  (19)  is  valid  for  the  time  period  (0  s  t  -  t£)  when  the  vapori¬ 
zation  process  follows  the  large  pai'ticle  model  (equation  (10A) ).  During 
the  time  period  when  the  vaporization  process  follows  the  small  particle 
model  (equation  (20))  ,  the  total  number  of  free  atoms  (neglecting  ioniza¬ 
tion)  is  the  sum  of  those  from  the  large-particle  vaporization  process 
which  occurred  earlier,  i.c.,  K[r03  -  rc3],  and  those  from  the  continuing 
small  particle  vaporization  described  by  equation  (15).  Note  that  it  is 
necessary  to  redo  fine  the  time  scale  used  in  that  portion  of  equation  (20) 
which  do.al.s  with  the  small  particle  vaporization  because  the  small  particle 
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vaporization  model  is  not  employed  until  after "the  completion  of  the  large 
particle  model,  i.e.,  t  >  t^.  Therefore,  the  time  becomes  t  -  t  in  equation  (20) 
when  the  free  atoms  are  being  formed;  however,  in  the  ionization  term, 
e  time  is  still  defined  from  the  beginning  of  the  initial  particle 

vaporization  because  ionization  decreases  the  total  free  atom  population 
at  all  times.  After  the  completion  of  vaporization,  t  >  t£,  equation  (21) 
will  be  valid  as  long  as  the  ionization  process  is  first-order  and  the  ion- 
electron  recombination  process  is  negligible.' 

EXPERIMENTAL 

Wl'lAAM/WW. 

To  study  the  vaporization  process,  emission  intensity-time  profiles 
were  obtained  using  the  droplet  generator  system  originally  described  by 
llieftje  and  Malmstndt:  [3].  In  this  technique,  droplets  are  individually 
and  reproduce bl y  injected  into  an  air-acetylene  flame  in  such  a  manner  that 
one  views  an  individual  droplet  and  its  products  without  interference  or 
influence  from  products  of  neighboring  droplets. 

The  air-acetylene  flame  was  supported  on  a  modified  Meker  burner 
[1]  and  was  operated  at  an  air  flow  of  15.0  L/min  and  an  acetylene  flow 
of  2.25  L/min.  The  velocity  of  the  flame  was  determined  to  be  11  m/s  [7]. 

A  quartz  sheath,  17.5  cm  above  the  burner  top,  was  used  to  stabilize  the 
fjame  [0].  The  plume  of  emitting  vapor  products  was  focused  without  magni¬ 
fication  onto  the  entrance  slit  (200  pm)  of  a  0.35  m  monochromator  (Model  EU  700, 
Heath  Co.,  Benton  Harbor,  MI),  which  was  tipped  on  end  so  the  entrance  slit 
i tit  e 1 1  •  m  ,  d  n  t  I *' t  o  lioi  i  .-.on i  a  I  . ions- sec l i on  of  the  flame  image.  The 
ini  i  i  on  .inii.il  w.o.  tit  ■  i  iv  1 1  .1  by  a  i:1'  ’S  phot  own  1 1  ipl  ier  tube  (llammamat  su  I'orp.  , 
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Middlesex,  N.J)  operated  at  -800  V  and  the  resulting  current  measured  by  a 
picoaminetcr  (414S,  Keithley  Instruments,  Inc.,  Cleveland,  OH).  Data  from 
the  picoaminetcr  were  collected  by  a  PDP  12/40  computer  (Digital  Equipment 
Corp. ,  Maynard,  MA) .  To  eliminate  the  effects  of  flame  temperature  variations, 
the  optical  system  was  fixed  in  position  to  monitor  the  emission  signal  at 
a  constant  height  in  the  flame  (19  cm  above  the  burner  top) .  To  obtain  the 
temporal  behavior  of  the  emission  intensity,  the  droplet  generator  was 
scanned  vertically  under  computer  control  by  means  of  a  stepper  motor  (SM-2A, 
Dcnco  Research,  Inc..,  Tucson,  A2)  and  a  lead  screw  connected  to  the  drop¬ 
let  generator  mount. 

Solutions  of  the  chlorides  were  prepared  from  dried  reagent  grade 
chemicals  according  to  standard  procedures  [9]. 

Droplet  diameters  were  measured  by  the  MgO  impression  technique  [10], 

The  initial  dried  particle  radius,  r  ,  was  calculated  indirectly  from  the 
initial  droplet  diameter,  solution  concentration  (64  pg/mL  in  all  cases),  and 
the  crystalline  solute  density  (NaCl,  2.169  g/mL;  KC1,  1.984  g/mL;  CsCl, 

3.988  g/mL) . 

RESULTS  AND  DISCUSSION 

'WW»flAAA.V\,VWVWVW\A. 

F-'igi  re  1  shows  the  experimental  emission  intensity-time  curve  obtained 
for  a  64  pg/mL  potassium  solution.  To  qualitatively  describe  this  curve  it 
is  convenient  to  divide  it  into  three  regions  as  shown  in  the  figure: 

Region  A)  At  the  beginning  ol  the  vaporisation  process,  the  emission  inten¬ 
sity  increases  rapidly  (vaporisation  rate  is  proportional  to  Ar2)  as  free 
aloms  and  other  species  are  vo  I  :<1  i  I  i  ned  from  the  surface  of  the  particle. 


....  > 
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Ionization  is  also  no nrr i nj;  in  this,  region  and  will  serve*  to  decrease  the 
emitting  atom  popu  I  at  ion .  Iteg  ion  H)  In  the  second  region,  t  he  vaporization 
process  continues  (the  vaporization  rate  is  now  proportional  to  Ar)  to  pro¬ 
duce  free  atoms  and  other  species  while  ionization  continues  to  remove  the 
emitting  atoms.  The  combination  of  the  slower  production  and  the  continual 
removal  of  the  free  atoms  causes  a  slowing  in  the  rate  of  growth  of  the 
emission  intensity.  Region  C)  At  the  completion  of  the  vaporization  process 
there  is  the  continued  removal  of  the  free  atoms  through  the  ionization  pro¬ 
cess.  This  results  in  a  decrease  in  the  emission  intensity. 

In  order  to  fit  equations  19-RI  to  an  experimental  curve,  it  is  nec¬ 
essary  to  obtain  values  for  K'  k  ,,  It.,  t  ,  r  ,  k  and  t...  These  values 

>  v l  1  c  c  v2  f 

are  calculated  from  five  data  point;;  selected  from 

the  experimental  curve  (cf.  ,  for  example,  Fig.  1).  The  first  data  point 
is  located  at  the  beginning  of  vaporization  (t  =  0,  1=0).  The  second  point 
is  chosen  on  the  rapidly  rising  portion  of  the  curve,  near  the  beginning  of 
vaporization.  A  third  point  is  selected  near  the  maximum  intensity.  The 
fourth  and  fifth  points  arc  chosen  a  short  time  beyond  the  maximum  intensity 
where  vaporization  is  presumed  to  be  complete.  Only  the  choice  of  data 
point  one  is  critical  (r  =  0,  1  =  0);  the  others  are  used  only  to  calcu¬ 
late  appi oxiinate  values  of  the  desired  parameters  and  establish  approximate 
time  zones  (large  particle  behavior,  small  particle  behavior,  time  after  com¬ 
pletion  of  vaporization)  which  are  necessary  to  define  the  regions  for 
these  calculations.  An  outline  of  the  method  is  described  below. 

After  vapoi i -at  ion  is  complete  (region  C) .  In  this  region,  where  only 
ionization  is  occur  ring,  k .  and  k* c.n  be  calculated  from  experimental  data 


point:.  -I  and  !>  using  equat  ion  (.11): 
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hi  -  [1/(1,  -  ts)]  ln( I5/U)  t5  >  t,  >  tf 


(22) 


(I,/ro3)e 


kit. 


t,  >  t. 


(23) 


Large  particle  time  period  (region  A) .  In  this  region,  kyl  can  be 
calculated  from  experimental  data  point  2  using  equation  (19)  and  the  values 
calculated  from  equations  (22)  and  (23). 
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0  <  t2  <  t 


(24) 


Small  particle  time  period  (J^gion  B) .  To  calculate  ky2 ,  a  value  for 
i*c  is  necessary.  The  critical  radius  is  conveniently  set  initially  at  zero 
(large  particle  model)  and  iteratively  increased  to  yield  the  best  calculated 
fit  to  the  experimental  data  points.  A  value  of  t  can  now  be  calculated 
from  equation  (17).  With  this  value  of  t£  and  experimental  data  point  3, 
it  is  possible  to  calculate  first  ky2  (equation  (25))  and  t^  (equation  (18)) 
using  equation  (20): 


v?. 
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(t3 
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tc  <  t3 


(25) 


Two  things  should  be  noted  concerning  the  method  of  calculating  the 
unknown  parameters  as  described  above:  (1)  the  experimental  data  point(s) 
used  to  calculate  a  specific  parameter  must  fall  within  the  time  zone  where 
that  calculation  is  valid;  (2)  once  a  value  has  been  calculated  for  each 
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of  the  curve- fit  parameters  amt  the  time  zones  established,  any  desired 
number  of  data  points  (we  normally  use  IS  for  the  entire  curve)  can  be 
assigned  to  the  appropriate  time  zones  and  the  corresponding  curve-fit 
parameter  can  be  calculated  and  averaged  for  all  data  points  that  fall 
within  that  time  zone. 

Figure  2  demonstrates  the  improvement  in  curve-fit  which  the  combined 
model  provides  over  either  the  large  or  small  particle  model  for  droplets 
of  a  64  gg/ml,  potassium  solution.  The  various  curve-fit  parameters  for  each 
of  the  three  models  based  on  I  ig.  2  are  shown  in  Table  1.  The  difference 
in  intensity  between  each  experimental  data  point  and  the  calculated  value 
is  squared  and  summed  for  all  of  the  experimental  data  points.  This  value 
is  reported  as  the  T.  (residual)'  in  Table  1.  The  sum  of  the  squares  of  the 
fsiduals  for  the  combined  model  indicates  dramatically  improved 
agreement  with  observations  over  Loth  the  large  and  small  particle  models. 

In  Figs.  3  and  4,  experimental  data  points  and  the  combined-model 
curve  are  shown  for  vapor i zat ton  of  a  particle  derived  from  solutions  of 
64  ug/mL  sodium  and  64  ug/ial.  cesium,  respectively.  The  combined-model  para¬ 
meters  for  these  emission  intensity  profiles  are  compiled  in  Tabic  2.  The 
accuracy  of  vaporization  rates  cited  in  Tables  1  and  2  is  limited  by  the  use 
of  approximate  solut  ions  to  di  1  Tena  nt  ial  equat  ions  (6)  and  (12).  In  a  com¬ 
parison  of  these  results  with  an  exact  numeric  solution  of  the  differential 
equations,  wo  have  found  an  error  in  calculated  vaporization  rates  up  to  10*. 
However,  the  simplicity  of  the  present  treatment  and  its  utility  for  calculating 
ionization  rates  renders  such  a  possible  error  acceptable  in  many  cases.  Work 
is  currently  underway  to  evaluate  numerically  t he  vaporization  rate  constants 


for  each  of  the  a (kali  met  a  I s  . 


Importantly,  the  error  associated  with  the 
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present  approximation  is  determinate,  and  will  be  greatest  for  metals  which 
exhibit  the  most  ionization  (e.g.,  Cs). 

Few  of  the  parameters  calculated  form  this  model  and  reported  in  Tables 
1  and  2  can  be  compared  with  literature  values  because,  for  the  most  part, 
such  values  do  not  exist.  The  values  of  the  ionization  rate  constants 
reported  here  are  in  rough  agreement  with  those  reported  in  other  flames 
[11].  The  values  of  the  particle  diameters  corresponding  to  the  critical 
radii  yielding  t lie  best  curve  fits  for  sodium,  potassium  and  cesium  arc 
0.88  pm,  0.86  pm,  and  0.64  pm,  respectively.  Using  appropriate  values  of 
the  diffusion  coefficients  for  the  elements  as  reported  by  SNELLEMAN  [12], 
values  of  the  mean  free  path  for  sodium,  potassium  and  cesium  are  calculated 
to  be  0.77  pm,  0.55  pm,  and  1.29  pm,  respectively  [4].  These  differences 
between  critical  diameters  and  mean  free  paths  could  be  due  r0  some  un¬ 
certainty  in  the  values  of  the  diffusion  coefficients  and/or  the  effect  of 
combining  the  large  and  small  particle  models  in  such  a  manner  that  the 
change  from  one  model  to  the  other  is  abrupt  in  our  treatment  whereas  in  a 
real  system  the  change  would  be  gradual . 


CONCLUSION 

In  this  study  a  combined  vaporization-ionization  model  for  single 
solute  particles  in  flames  is  described.  This  model  provides  a  new  method 
for  the  determination  of  the  rates  of  ionization  and  particle  vaporization. 
For  the  first  time,  approximate  vaporization  rate  constants  for  three  alkali 
chlorides  arc  reported.  Although  the  model  has  only  been  applied  to  a 
selected  number  ,»f  elements,  it  is  believed  to  be  of  general  utility  for 
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systems  which  exhibit  first  onler  ionization.  The  model  can,  however,  cause 
large  systematic  errors  in  calculated  vaporization  rates  when  the  ioniza¬ 
tion  rates  are  large.  Work  with  exact  numeric  solutions  of  the  differential 
equations  are  being  examined  in  more  detail. 
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Calculated  based  on  solution  concentration,  an  assumed  density  of  1.984  g/mL,  and  an  initial  droplet 
diameter  of  27. S  urn. 
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figure  1. 


1-xperimental  data  points  and  the  combined  model  curve-fit 
for  vaporization  of  particles  obtained  from  a  64  pg/mL  cesium 
solution.  Initial  droplet  diameter  is  27.5  pm.  □  are  the 
experimental  data  points. 
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